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Preface 


1  For 


Our  work  during  the  period  covered  by  this  report  was 


devoted  to  the  properties  of  materials  using  methods  of  statis-  a  ^ 

'  ion _ 

tical  mechanics.  - - - - - 

The  main  focus  of  our  research  was  in  the  following  ■ 

areas: 

,  Dist  /  '"sds  an5/0r 

t.i  Kinetics  of  quenched  alloys.  peoiai 


tl'/.y  Phase  diagrams  of  alloys.  IfcLLl 

III.  Equilibrium  and  kinetic  properties  of  polymers, 

IV*  •  Properties  of  simple  fluids  and  electrolytes. 

Our  work  involved  both  theoretical  studies  and  computer 
simulations  of  model  systems.  These  simulations,  which  were 
among  the  largest  of  their  kind,  greatly  contributed  to  our 
understanding  of  the  microscopic  processes  involved  in  various 
macroscopic  phenomena.^ 

I.  ALLOYS 

1.  Kinetics  of  Phase  Transformations  in  Quenched  Alloys. 

During  the  whole  period  covered  by  this  report  we  were 
engaged  in  a  multifaceted  study  of  the  time  evolution  of  quenched 
binary  (A-B)  alloys.  At  sufficiently  high  temperatures  all 
alloys  exist  in  a  single  phase  disordered  state.  At  low  tempera¬ 
tures,  they  generally  exist  in  ordered  states.  The  nature  of  the 
ordered  states  depends  on  whether  it  is  preferable  (energetically) 
for  the  particles  of  each  species  to  be  surrounded  by  their  own 
kind  or  by  the  different  species:  examples  of  such  systems  are 

ZnAl  and  CuPd,  respectively.  ($fc, 
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ferromagnetic  and  anti-ferromagnetic  states  in  magnetic  systems. 

The  question  of  interest  is  how  does  an  alloy  which  is 
suddenly  cooled  (quenched) ,  so  that  the  atoms  have  no  time  to 
migrate  during  the  quench,  evolve  in  time  (following  the  quench) 
from  the  homogeneous  disordered  state  it  finds  itself  in  to  the 
equilibrium  phase  appropriate  for  it  at  the  low  temperature; 
phase-segregated  or  ordered.  This  is  a  problem  of  practical 
importance  to  the  metallurgist  since  properties  of  many  alloys 
such  as  their  strength,  ductility,  etc.  depend  on  their  struc¬ 
ture.  This  in  turn  is  greatly  influenced  by  the  heat  treatment 
the  material  undergoes.  An  tinderstanding  of  the  kinetics  of 
such  processes  should  thus  aid  in  the  development  of  materials 
with  desirable  properties, 
a)  Segregating  Alloys 

We  studied  the  structure  functions  S(k,t)  obtained 
from  computer  simulations  of  the  time  evolution  of  a  segregat¬ 
ing  binary  alloy  following  quenching  into  the  immiscibility  gap. 

Our  most  important  finding  was  that  S(k,t)  has  a  simple  scaling 

s 

behavior,  S(k,t)^K  (t)F(k/K(t))  at  late  times.  The  shape  of  the 
function  F(x)  depends  somewhat  on  the  part  of  the  coexistence 
region  into  which  the  quench  is  made.  Comparison  with  experi¬ 
ments  on  quenched  alloys  and  also  binary  liquids  is  quite  satis¬ 
factory.  In  fact  most  experimentalists  now  use  our  results  as  a 
standard  of  comparison  for  analyzing  their  data:  see  publication 
[25]  for  the  most  complete  description  of  our  results. 
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We  also  studied  the  time  evolution  of  the  "grains" , 
i.e.  clumps  of  precipitate  of  the  minority  phase,  following 
quenching.  Some  features  of  the  size  distribution  of  the 
grains  were  found  to  satisfy  theoretical  predictions  others 
did  not.  The  relation  between  grain  size  distribution  and 
structure  function  were  investigated  and  some  aspects  clari¬ 
fied.  Further  theoretical  and  computer  experiments  are  needed 
to  obtain  a  full  understanding  of  the  coarsening  process, 
b)  Ordering  Alloys 

Alloys  which  form  ordered  superlattices  at  low 
temperatures  are  described  quantitatively  by  a  long-range- 
order  parameter  n,e<1  which  vanishes  for  T>TC  and  is  non- 
zero  for  T<TC« 

In  our  work  we  carried  out  the  first  quantitative 
computer  simulations  of  the  ordering  process  in  a  model  binary 
alloy  following  quenching  from  a  disordered  (usually  molten) 
state  to  a  temperature  T<TC. 

We  analyzed  our  simulation  data  on  the  basis  of  the 
movement  of  domain  walls  separating  regions  (domains)  whose 
ordering  is  different.  The  motion  of  the  domain  walls  is  such 
as  to  reduce  their  curvature.  The  theory  predicts  that  the 
surface  area  per  unit  volume  will  decrease  in  proportion  to 
The  results  of  our  computer  simulations  gave  strong 
support  to  the  theory.  We  found  the  predicted  t law  dur¬ 
ing  early  periods  following  the  quench,  when  the  characteristic 
length  £  of  the  domains  is  small  compared  to  the  size  of  our 
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system.  Also  during  this  period  the  "staggered”  structure 
function  S(k,t)  is  observed  to  scale  in  the  way  the  theory 
suggests.  When  5  becomes  comparable  to  the  size  of  the 
system,  however,  finite-size  effects  become  important  and 
there  are  large  differences  in  the  behavior  of  different 
samples  in  our  simulation.  Nevertheless,  the  late  stages  of 
the  process  can  also  be  interpreted  in  terms  of  domain  walls. 

The  physical  situation  here  is  quite  different  from 
that  in  a  segregating  alloy  in  which  the  order  parameter  is 
conserved.  While  there  too,  we  found  scaling  behavior,  the 
physical  interpretation  (and  the  scaling  exponent)  are  quite 
different. 

II.  PHASE  DIAGRAMS  OF  BINARY  ALLOYS 

We  carried  out  extensive  Monte  Carlo  computations  of 
the  coherent  phase  diagram  (in  the  temperature-composition 
plane)  of  ordering  binary  alloys  on  a  face-centered  cubic 
lattice.  Long-  and  short-range  order  parameters  as  well  as 
ordering  energies  were  obtained.  We  considered  systems  with 
both  nearest  and  next  nearest-neighbor  interactions.  Our  re¬ 
sults  were  compared  with  various  theoretical  predictions,  parti¬ 
cularly  with  the  cluster- variation  method  which  was  found  to  be 
the  most  reliable.  It  gives  fairly  accurate  results  at  stoich¬ 
iometric  composition.  However,  it  too  was  found  to  fail  in  the 
more  general  case  predicting  a  wrong  topology  of  the  phase 
diagram  in  the  case  of  nearest-neighbor  interaction.  Our  work 
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still  stands  now  as  the  most  reliable  "experimental"  phase 
diagram  for  these  systems. 

III.  POLYMERS 

When  a  solution  of  macromolecules  is  sufficiently  dilute, 
the  interactions  between  polymer  chains  can  be  neglected.  A 
good  model  for  this  system  then  consists  of  one  chain  in  a  sea 
of  solvent . 

In  our  work  we  first  investigated  theoretically  and  via 
computer  simulations  the  statics  and  dynamics  of  a  polymer  chain 
in  a  "good"  solvent  where  only  excluded  volume  effects  are  impor¬ 
tant.  We  found  that  the  distribution  function  for  two  beads 
scales  as  the  number  of  beads  N  to  the  power  v  and  that  edge 
effects  are  small.  The  dynamical  correlation  functions,  such  as 
that  of  the  end-to-end  vector,  scale  as  N2v+1  with  v  *  0.6. 

The  results  of  a  dynamical  lattice  polymer  model  were  shown  to  be 
consistent  with  our  results  if  one  adjusts  the  time  scales  in 
such  a  way  that  the  center  of  mass  diffuses  at  the  same  rate  in 
the  two  models.  The  relaxation  of  the  stress  tensor  was  found  to 
be  quite  similar  to  that  of  the  Rouse  model.  Finally,  it  was 
shown  that  edge  effects  are  much  more  pronounced  in  the  diffusive 
motion  of  the  individual  beads,  there  being  a  skin  comprising 
about  40%  of  the  total  polymer,  where  bead  motion  is  relatively 
quicker. 

Our  results  confirm  the  idea  that  for  long  chains  the 
degree  of  expansion  due  to  excluded  volume  depends  on  the  inter- 


action  range  o,  on  the  number  of  links  N,  and  on  the  link  size  b 
through  a  single  variable  z  «  (a/ b)  /“JT. 

We  then  studied  polymer  chains  in  a  "poor  solvent”  with 
both  repulsive  and  attractive  interactions.  Varying  the  tempera¬ 
ture  represents  a  polymer  chain  in  solvents  with  differing 
degrees  of  solubility.  We  found  that : 

(a)  At  high  temperatures,  corresponding  to  good  solvents,  the 
dependence  of  the  chain  dimension  R  on  the  number  of  units 

N  follows  the  same  power  law  as  does  a  chain  with  purely  re¬ 
pulsive  interaction;  R^  *  N^u,  u  *  0.6. 

(b)  There  is  a  temperature  If  in  the  neighborhood  of  Flory's 
0-temperature  (at  which  the  second  virial  coefficient  of 
the  interaction  vanishes)  where  there  is  a  linear  relation 

A, 

between  R^  and  N:  0-0  depends  on  the  "stiffness"  of  the 
chain . 

(c)  At  lower  temperatures  the  chain  manifests  a  collapse  which 
becomes  more  pronounced  with  increasing  N. 

There  is  overall  qualitative  agreement  between  our  results,  experi¬ 
ment,  and  a  generalized  form  of  Flory's  theory  although  the  behavior 
in  cases  b)  and  c)  appears  to  depend  on  more  details  of  the  inter¬ 
action  than  are  taken  into  account  in  the  latter. 

We  also  carried  out  computer  simulations  on  a  polymer  chain 
with  excluded  volume  interactions  in  the  presence  of  an  external 
stretching  force.  For  weak  and  moderate  forces  the  response  is 
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linear  while  for  strong  forces  the  behavior  is  nonlinear — 
consistent  with  the  non-Gaussian  nature  of  the  end-to-end 
vector  R  distribution  for  large  R.  In  the  vicinity  of  the 
e-temperatures  the  onset  of  nonlinearity  occurs  at  larger 
forces. 

Polymer  chains  with  excluded  volume  interactions  which 
are  confined  within  a  slab  of  width  D  were  also  investigated. 

We  studied,  in  particular,  R2(N,  D) ,  the  mean  square  of  the 
end-to-end  distance  of  a  chain  with  N  links.  For  large  N  and 
fixed  D,  R2(N,D)  r|(N)  ,  the  squared  end-to-end  distance  of  a 

o  p  2\) 

chain  constrained  to  a  plane.  We  found  R^CN):  R  (N,0)  «  N 
with  2v  »  3/2  in  agreement  with  the  prediction  of  Flory.  Letting 
r|(N)  =  R2(N,®),  the  end-to-end  distance  for  an  unconstrained 
three  dimensional  chain,  we  examined  the  cross-over  scaling  of 
R2(N,D)/R2(N)  as  a  function  f(x)  of  x  s  D/R3(N).  For  x  <  0.45, 
f(x)  *  x"1/2,  in  agreement  with  predictions  of  Daoud  and  de  Gennes 
IV.  SIMPLE  FLUIDS  AND  ELECTROLYTES 

Our  interest  in  the  equation  of  state  and  correlation 
functions  of  fluids  dates  back  a  long  time  and  during  the  period 
covered  by  this  report  we  continued  some  work  in  this  area.  In 
particular  we  investigated  the  structure  of  simple  liquids  as 
well  as  that  of  electrolytes  in  the  vicinity  of  a  wall.  (For 
electrolytes  the  wall  can  represent  a  solid  electrode  and  may 
therefore  carry  an  electric  charge).  Our  work  was  two  pronged: 
a)  we  looked  for  and  found  some  exact  sum  rules  which  gave  exact 
relations  for  suitably  weighted  integrals  over  the  correlation 


functions,  b)  we  used  various  approximation  schemes,  e.g.  Percus- 
Yevick  integral  equation,  Mean  Spherical  Approximation  and  its 
generalizations,  to  compute  the  full  density  profile.  The  two 
approaches  complemented  each  other  in  that  we  could  use  the  sum 
rules  for  fixing  certain  parameters  in  the  approximations. 

Our  results  during  this  period  included:  a)  an  exact 
general  expression  for  the  contact  value,  p(0),  of  the  density 
profile  of  the  fluid  near  a  wall  which  could  be  charged,  b)  an 
exact  result  for  the  high  field  limit  of  the  contact  value,  q(0), 
of  the  charge  profile,  c)  full  density  profiles  in  good  agreement 
with  machine  calculations. 
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